Local Spin Anisotropy Effects upon the Magnetization and Specific Heat of Dimer 

Single Molecule Magnets 
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We present an exactly solvable model of equal spin Si dimer single molecule magnets. The 
spins within each dimer interact via the Heisenberg and the most general quadratic global and local 
(single- ion) anisotropic spin interactions, and with the magnetic induction B. For antiferromagnetic 
couplings and si > 1/2, the low temperature T magnetization M(B) exhibits 2si steps of universal 
height and midpoint slope, the sth step of which occurs at the non-universal level-crossing magnetic 
induction Bg 31 (0, cf>), where 0,<j> define the direction of B. The specific heat Cv exhibits zeroes 
as T — > at these Bf si (6, cj>) values, which are equally surrounded by universal peak pairs as 
T — > 0. The non-universal Bg tS1 {0, </>) values lead to a rich variety of magnetization plateau behavior, 
the structure and anisotropy of which depend upon the various global and local anisotropic spin 
interaction energies. We solve the model exactly for si = 1/2, 1, and 5/2, and present M(B) 
and Cv(B) curves at low T for these cases. For weakly anisotropic dimers, rather simple analytic 
formulas for M(B) and Cv(B) at arbitrary si accurately fit the exact solutions at sufficiently low T 
or large B. An expression for B^ S1 (6, 4>) accurate to second order in the four independent anisotropy 
energies is derived. Our results are discussed with regard to existing experiments on si = 5/2 Fe2 
dimers, suggesting further experiments on single crystals of these and some si = 9/2 [Mn^ dimers 
are warranted. 

PACS numbers: 05.20.-y, 75.10.Hk, 75.75.+a, 05.45.-a 



INTRODUCTION 

Single molecule magnets (SMM's) have been under in- 
tense study recently, due to their potential uses in mag- 
netic storage and quantum computing. [1-3] The mate- 
rials consist of insulating crystalline arrays of identical 
SMM's 1-3 nm in size, each containing two or more 
magnetic ions. Since the magnetic ions in each SMM 
are surrounded by non-magnetic ligands, the intermolec- 
ular magnetic interactions are usually negligible. Al- 
though the most commonly studied SMM's are the high- 
spin Mni2 and Fes, [1-5] such SMM's contain a variety 
of ferromagnetic (FM) and antiferromagnetic (AFM) in- 
tramolecular interactions, rendering unique fits to a va- 
riety of experiments difficult. [6] 

In addition, there have been many studies of AFM 
Fe„ ring compounds, where n = 6,8,10,12, etc. [7-10] 
In these studies, analyses of inelastic neutron diffraction 
data and the magnetic induction B dependence of the 
low-tempcrature T specific heat and magnetization steps 
were made, using the isotropic Heisenberg near-neighbor 
exchange interaction, the Zeeman interaction, and var- 
ious near- neighbor spin anisotropy interactions. [7-10] 
However, the rings were so complicated that analyses of 
the data using those simple models were inaccessible to 
present day computers. [8, 9] Thus, those authors used ei- 
ther simulations or phcnomcnological fits to a first-order 
perturbation expansion with different spin anisotropy 
values for each global ring spin value. [8-10] 

Here we focus on the simpler cases of equal spin s\ = s% 
magnetic dimers, for which the full spin anisotropy effects 



can be evaluated analytically, investigated in detail nu- 
merically, and compared with experiment. AFM dimers 
with si = 1/2, 3/2, [11-14] and various forms of Fe2 with 
si = 5/2 were studied recently. [15-20] Several Fe2 dimers 
and effective si = 9/2 dimers of the type [Mn4]2,[21, 22] 
have magnetic interactions weak enough that their ef- 
fects can be probed at T « IK with presently available 
B. A comparison of our results with magnetization M 
versus B step data on a Fe2 dimer strongly suggests a 
substantial presence of local spin anisotropy. [17] 

The paper is organized as follows. In Section II, we 
present the model in the crystal representation and given 
exact formulas for the matrix elements. The general ther- 
modynamics are presented in Section II, along with the 
universal behavior of the M(B) and Cv{B) behavior as- 
sociated with the energy level crossing. In Section IV, we 
solve the model exactly for s\ = 1/2, giving analytic ex- 
pressions for the magnetization and specific heat. In Sec- 
tion V, we discuss the exact solution for si = 1, present 
the equations from which the eigenvalues are readily ob- 
tained, and give numerical examples of the low-T mag- 
netization and specific heat curves. In Section VI, we 
present numerical examples of the low-T s± — 5/2 mag- 
netization and specific heat curves. In section VII, we 
rotate to the induction representation, and give the eigen- 
states to first order in the anisotropy energies. These are 
used to obtain asymptotic expressions for the magnetiza- 
tion and specific heat for arbitrary si that are highly ac- 
curate at sufficiently low T and/or large B. In addition, 
analytic formulas for the magnetic inductions at which 
the level crossings occur are provided, accurate to sec- 
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ond order in the anisotropy energies. Finally, in Section 
VIII we discuss our results with regard to experiments 
on Fe 2 dimers. 



THE MODEL IN THE CRYSTAL 
REPRESENTATION 

We represent the s\ = s 2 dimer quantum states, \ip™} 
in terms of the global (total) spin and magnetic quantum 
numbers s and m, where S = Si + S 2 and S z = S ■ z 
satisfy S 2 \4>™) = s(s + l)\ip™) and S z \ip™) = m\ip™), 
where s = 0, 1, . . . , 2s\, m = — s, . . . , s, and we set K = 1. 
We also have S±\ip™) = A± m |?A™ ±1 ), where S± = S X ± 
iS y and 



= sj (s - m)(s + m + 1). 



(1) 



For an arbitrary B, we assume the Hamiltonian has the 
form H — Hq + H a + Hb + Hd + H e , where 



Ho = -JS 2 /2- 1 S-B 



(2) 



contains the Heisenberg exchange and Zeeman interac- 
tions, the gyromagnetic ratio 7 = g^B, where g ~2 and 
Hb is the Bohr magneton. The global axial and azimuthal 
anisotropy terms 



Hb 



JhS 2 



and 



Hd — —Jd(Sl - S 2 ), 



(3) 



(4) 



respectively, only involve components of S, but have 
been the main anisotropy terms discussed in the SMM 
literature, [5, 21] so we have included them for compari- 
son. Such terms have been commonly studied as an ef- 
fective Hamiltonian for a singlet orbital ground state, in 
which the tensor global spin interaction with a fixed spin 
quantum number s has the form S ■ A • S, resulting in the 
principal axes x,y and z. [23] For a dimer, we take A to 
be diagonal in the orientation pictured in Fig. 1. Then 
J = -(A xx + A yy )/2, J b = -A zz + (A xx + A yy )/2, and 
Jd = (A yy - A xx )/2. Taking \ J d /J\ < 1 and \ J b /J\ < 1 
still leaves Jd/ Jb unrestricted. The single-ion axial and 
azimuthal anisotropy terms, 



H a — J a Sf. 



and 



H e 



i=l 
2 

t=l 



(5) 
(6) 



respectively, arise from spin-orbit interactions of the local 
crystal field with the individual spins. These terms have 
usually been neglected in the SMM literature, but have 



been studied with regard to complexes containing a single 
magnetic ion, such as Ni +2 ,[23] and with regard to clus- 
ters of larger numbers of identical magnetic ions. [23, 24] 
The local axially and azimuthally anisotropic exchange 
interactions 



satisfy 



Hf — —JfSi z S2 Z , 

H c = —J c {Sl x S2 x — S\ y S2y), 



2Hf/Jf = Hb/Jb-Ha/Ja, 
2H C /J C = Hd/Jd-He/Je, 



(7) 
(8) 



(9) 
(10) 



so we need only include either either H a or Hf and H e 
or H c , respectively. [23] That is, if we stick to the Hamil- 
tonian H, we may incorporate the effects of Hf and 
H c by letting J b — > J b + J//2, J a — > J a - J//2, and 
Jd — ► Jd + Jc/2, J e — > Je — Jc/2, respectively. Since H a 
and H e describe the axial and azimuthal anisotropy each 
single ion attains from its surrounding environment, they 
are the physically relevant local anisotropy interactions. 

For the case of Fe2,[15] a constituent of the high-spin 
SMM Fe 8 and the AFM Fe„ rings, [4, 8, 9] the exchange 
between the Fe +3 si — 5/2 spins occurs via two oxygen 
ions, and these four ions essentially lie in the same (xz) 
plane. [15, 18] We set the z axis parallel to the dimer axis, 
as pictured in Fig. 1. Since the quantization axis is along 
the dimer axis, which is fixed in a crystal, we denote this 
representation as the crystal representation. 

We generally expect each of the Jj for j = a, b, d, e to 
satisfy \Ji/J\ <SC 1, but there are not generally any other 
restrictions upon the various magnitudes of the Jj . Since 
all dimers known to date have predominantly AFM cou- 
plings ( J < 0), and also because their magnetizations and 
specific heats are particularly interesting, we shall only 
consider AFM dimers. In addition, since no studies on 
unequal-spin dimers have been reported to our knowl- 
edge, we shall only treat the equal-spin s\ — s 2 case. 
We note that for equal spin dimers, the group symmetry 
of the dimer environment is C 2v , so that Dzaloshinskii- 
Moriya interactions do not arise. [23] Hence, our Hamil- 
tonian H is the most general quadratic anisotropic spin 
Hamiltonian of an equal-spin dimer. 

For B = B{sin 9 cos </>, sin 6 sin 0, cos 6) , we have 

WolO = E?\W)+5E ]T e-^AHV> s ro+CT >, 



cr=±l 



and 



HdW) = 
where 



V F 
2 ^ 

<7=±1 



(11) 

(12) 



(13) 
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FIG. 1: Sketch of an Fe2 dimer, with two bridging O -2 ions 
(O). Ligands (not pictured) are attached to the Fe +3 ions 
(Fe). The arrows signify spins. 



-Js(s + l)/2 - mbcos9, 

-h 
2 

■yB. 



-b sin 



E T 

SE 
b 

and 



Ff = A x s Al +x . 



(14) 
(15) 
(16) 

(17) 



7i a and Ti. c contain the individual spin operators Si Z 
and Si± for i = 1,2. After some Clebsch-Gordon alge- 
bra involving the Wigner-Eckart theorem, for arbitrary 

(sx,s,m), 

-cr^Timf)), (18) 

777 I / 

«r) = ^w)-^-i)\d:: s1 w-i) 



where 



and 



CZi = ^{s-m){s-m- rjr) s , Sl , 
D™ S1 = V(s 2 -m^ Sl , 



Vs, Sl = V[(2si + l) 2 -s 2 ]/(4 S 2 -l). 



(19) 

(20) 
(21) 

(22) 



For s = 0, we require m = 0, for which Cq = -Do, Sl = 0. 
We then find 



U e Ws 



G'' 



Je 

4 



E *C'f>. 

'=±1 

E (^Tj^; 
E ^T/iv> 



s+2a' 



m+2<x\ 



= ±1 



m-^-2(T , 
s+2cr' . 



(23) 



(24) 



a' = ±l 
s(s + 1)-1 



and 



UlTl.cr 
H s,si 

K x,a' 

8, Si 

T X 

^8,8! 



+[m J + l-2a(a + l)]a a , ai , 

— n m n m 

— iy s+(<j' + l)/2,s 1 i/ s+(3<T' + l)/2 : si' 

_ r ,-<j'a;-(l+ ( T')/2 rr -(j'a;-(3cr' + l)/2 

— Ly s+(cr' + l)/2,s 1 U s+(3<r' + l)/2,si ' 

= 2Ffa.,. s1 , 



3s(s + l)-4si(si + l)-3 



(2s-l)(2s + 3) 



(25) 
(26) 

(27) 
(28) 

(29) 



We note that a SySl = l — — and that Eq. 

(29) holds for s > 0. Equations (18) and (19) allow 
for an exact solution to the most general Hamiltonian of 
arbitrary order in the individual spin operators. In all 
previous treatments of more complicated spin systems 
with similar anisotropic interactions, it was only possi- 
ble to obtain numerical solutions, and therefore the full 
anisotropy of the magnetization and specific heat was 
not calculated. [23, 24] The operations of 7i , Ti-b and 
Tid satisfy the selection rules As = 0, Am = 0, ±1,±2. 
The local anisotropy interactions Ti. a and TL e allow tran- 
sitions satisfying As = 0,±2, Am = 0, and As = 0,±2, 
Am = ±2, respectively, so in the presence of either of 
these interactions, s is no longer a good quantum num- 
ber, unless si = 1/2. 



GENERAL THERMODYNAMICS 

In order to obtain the thermodynamic properties, we 
first calculate the canonical partition function, Z = 
Trexp(— f3H). Since TL is not diagonal in the (s,m) rep- 
resentation, we must construct the wave function from 
all possible spin states. We then write 



Z = Tr^Je-^l*^), 



(30) 



where is constructed from the {|^™)} basis as 

(31) 

where j3 = l/(fcsT) and fee is Boltzmann's constant. 
To evaluate the trace, it is useful to diagonalize the 
(y si \H\y si ) matrix. To do so, we let |* Sl ) = U\$ Sl ), 
where 



(32) 



is constructed from the new orthonormal basis {|^n)}, 
and U is a unitary matrix of rank n Sl = (2si + 
l) 2 . Choosing U to diagonalize H., we generally obtain 
r H\<t> n ) = £n\4>n} and the partition function for a SMM 
dimer, 



Z = ^ exp(-/3e r , 



(33) 



4 



The specific heat Cv — k B 2 d 2 \nZ/df3 2 is then easily 
found at all T, B, 



C v = 



k B (3 2 
Z 2 



zte 2 n e-^-(±e n e-^)' 



71=1 



n=l 



.(34) 



It is easily seen from Eq. (34) that at the induction 
B 1 g Si corresponding to the sth level crossing in which 
e s = e s _i, Cy — > as T — > 0. The magnetization 



M = -J2 V B (en)exp(-/3e n ), 



(35) 



n=l 



requires Vs(e„) for each £?. As T — > 0, at most two 
eigenstates are relevant. For most B values, only one 
e n is important. But near the sth level-crossing induc- 
tion B l s c Si (9,<j)) at which e s = e s _i, two eigenstates are 
relevant. We then have 



0. 



C v 



sJ; ai (M)± 



2c 

7^ 



peak 



(36) 
(37) 



fe B T<|J| 

w 0.439229, (38) 



and 



M 



7<£B 



07 
4 ' 



(39) 
(40) 



where s = 1, . . . , 2si and c w 1.19967864 is the solution 
to tanhc = 1/c. The easiest way to obtain Eq. (39) is to 
first rotate the crystal so the quantization axis is along B, 
as discussed in Appendix B. We note that B ] f si (6, <j>) de- 
pends upon s,si, and the direction of B when anisotropic 
interactions are present. Hence, the heights and midpoint 
slopes of the 2si M(B) steps are universal, but the step 
positions and hence their plateaus are not. Correspond- 
ingly, the 2sj positions of the Cv(B) zeroes are non- 
universally spaced, but each zero is equally surrounded 
by two peaks of equal, universal height, the universal 
spacing between which is oc T as T — > 0. Hence, the 
non-universal level-crossing inductions B\ c Si (9, <f>) fully 
determine the low-T thermodynamics of AFM dimers. In 
the next three sections, we consider the special cases of 
si = 1/2, 1 and 5/2. Then, in Section VII and Appendix 
D, we present out general expression for B l f Si (6, </>) ac- 
curate to second order in each of the Jj. We remark 
that a double peak in the low-T Cv(B) curve has been 
seen experimentally in a much more complicated Fe6 ring 
compound, and was attributed to level crossing. [8] 



ANALYTIC RESULTS FOR SPIN 1/2 

Plots of Cv/ks and M/j versus jB/\J\ for the 
isotropic spin 1/2 dimer were given previously. [14] For 
si = 1/2 with an arbitrary B and Jj for j = a,b,d,e, 
the rank 4 Hamiltonian matrix is block diagonal, since 
s = 0, 1 is a good quantum number. The eigenvalues are 
given by 



Ja 

2 ' 

Ja 

2 



J + A„, 



(41) 

n = 2,3,4, (42) 



where 



= \l + 2\ 2 l J b -\ n [J 2 - J 2 + b 2 } 

-b 2 sin 2 9[J h - J d cos(20)]. (43) 

For the special cases B\\i for i = x,y, z, the X l n satisfy 



\ x -y 



0,-J b ±F z , 



-2 J 



-7 + F 



where 



b 2 + J?, 



F 



Jx.y,z — (*^6 ^ Jd)/^-) Jdf 



(44) 
(45) 



(46) 
(47) 



respectively, and where J x (J y ) corresponds to the upper 
(lower) sign. 

The magnetization for B\\i is given by 



Mi = 



^BsmHpFi) 



(48) 
(49) 



FiVi 

V, = cosh(/3F 4 ) + Ai/2, 
A x . v = exp(-/3J Xj!/ )[exp(-/3J)+exp(2/3J Xi j / )],(50) 
A z = cxp(-/3J b )[exp(-/3J) + l], (51) 



respectively. When the interactions are written in terms 
of the less physical J/, J^, and J c , then Jt — > J5 + Jf 
and Jd — > J a + J c /2, so that for s\ = 1/2, Jf and J c 
merely renormalize J^ and J^. Neither of the single-ion 
spin anisotropy terms H a and H e affect the thermody- 
namics for si = 1/2. We note that M y (Jd) = M x (—Jd) 
for each B, as expected from Hd = — Jd{S 2 — S 2 ). 

From Eqs. (41), (42), (44), and (45), the single level 
crossing induction (s = 1) for am si = 1/2 dimer with 
B\\i occurs at 



7^1 



/2 



\J J 2 + J(Jb ± Jd) , B\\x,y 
V(J + J b ) 2 - Jl B\\z 



(52) 



provided that J + J x . y < and J + Jf, < 0, respectively. 

To distinguish the different effects of the global 
anisotropy interactions J b and Jd that affect the mag- 
netization of si = 1/2 dimers, in Figs. 2 and 3, we 



■5 




g^ B B/|J| 



FIG. 2: Plots of M/7 versus jB/\J\ at k B T/\J\ = 0.03 for 
the AFM spin 1/2 dimer with J b = 0.1J, J d = 0, with B\\z 
(solid), B\\x (dashed), along with the isotropic case J b — 
J d = (dotted). 




U1 =-> — : 1 1 1 

0.75 1.00 1.25 



9H B B/|J| 



FIG. 3: Plots of M/7 versus jB/\J\ at k B T/\J\ = 0.03 for 
the AFM spin 1/2 dimer with J d = 0.1J, J b = 0, with B\\z 
(solid), B\\x (dashed), along with the isotropic case Jb = 
J d = Q (dotted). 



have respectively plotted the low-T M/j versus jB/\J\ 
with J b = 0.1J, J d = and J d = 0.1J, J b = for B\\\z 
(solid) and B\\x (dashed), along with the isotropic case 
■h = Jd = (dotted). From Fig. 2, J b < and J d = 
causes a greater shift to higher B values at the magneti- 
zation step with B\\z than for B\\x, consistent with Eq. 
(52). In addition, J d finite with J b = has a very differ- 
ent effect upon the anisotropy of the magnetization step, 
as shown in Fig. 3. Although for B\\z, B at the step is 
slightly reduced from its isotropic interaction value, for 
B\\x, the magnetization step occurs at a larger B. These 
results are consistent with Eq. (52). The midpoint slopes 
are universal, in accordance with Eq. (??slope). 
The specific heat of an si = 1/2 dimer with B\\i is 



a 



Vi 



(53) 



where the X>j are given by Eq. (49), and the A/j are 
given in Appendix A. Plots at low T of Cy/fc^ versus 
7-B/l J\ for S\ = 1/2 dimers with the corresponding global 
anisotropics J b = 0.1J and J d = 0.1J are shown in Figs. 
4 and 5, respectively. We note the universal curve shapes, 
but non- universal level-crossing positions, in quantitative 
agreement with Eqs. (36)- (38). In Fig. 4, the positions 
of the maxima and the central minimum in Cv track 
that of the magnetization step in Fig. 2 with the same 
parameters. With J d = 0.1J, the behaviors in Cv and 
M for B\\x are also very similar. However, there is a 
slight difference in the behaviors for B\\z. Note that M 
(Fig. 3) shows a slight reduction for B\ \z in the induction 
required for the step, whereas Cv (Fig. 5) shows a slight 
increase in the positions of the peaks. This detail only 
appears when J d ^ 0, for which the effective temperature 
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FIG. 4: Plot of Cv/k B versus jB/\J\ for the AFM spin 1/2 
dimer with J b = 0.1J, J d = at k B T/\J\ = 0.03 with the 
same curve notation as in Fig. 2. 



is slightly higher than for J d = 0. 



ANALYTIC AND NUMERICAL RESULTS FOR 
SPIN 1 

For dimers with si = 1, the allowed s values are s = 
0,1,2. The three s = 1 states are decoupled from the six 
remaining s — 0, 2 states. They satisfy a cubic equation 
given in Appendix A. For B\\i with i — x, y, z, this cubic 
equation simplifies to a linear and a quadratic equation, 
as for the s — 1 eigenstates of si — 1/2 dimers. 

The remaining six eigenstates corresponding nominally 
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FIG. 5: Plot of Cv/k B versus -yB/\J\ for the AFM spin 1/2 
dimer with J d = 0.1J, J b = at k B T/\J\ = 0.03 with the 
same curve notation as in Fig. 3. 



to s = 0, 2 are in general all mixed. The matrix leading to 
the hcxatic equation from which the six eigenvalues can 
be obtained is given in Appendix A. That is sufficient 
to evaluate the eigenvalues for si = 1 using symbolic 
manipulation software. When combined with the three 
s = 1 eigenvalues, one can then use Eqs. (34) and (35) 
to obtain the resulting exact magnetization and specific 
heat at an arbitrary B. The combined nine eigenvalues 
depend upon all four anisotropy parameters Jj for j = 
a, b, d, e. 

To the extent that the eigenvalues can be obtained 
from the solutions to either linear or quadratic equations, 
the expressions for the B 1 ^ are simple, and are given in 
Appendix A. In the global anisotropy case J a — J e — 0, 
the first level crossing induction B l £ A with B\\i is iden- 
tical to -Bj ^, the level crossing with s\ = 1/2 given by 
Eq. (52). For comparison, we expand the first and sec- 
ond level crossing inductions to first order in each of the 
Jj for j — a, b, d, e for B\\i where i = x,y, z, 



7<S 

lc(l) 



J + 


Ja 

3 


Jb, 




(54) 


J - 


Ja 

6 


Jb l (T 

T T 2 {Jd 




(55) 


2J- 


Ja- 


- 3Jfc, 




(56) 


2J4 


Ja 

2 


Jb 5J d 
~ ~2 T ~2~ 


3J e 


(57) 



tB 1c(1) 



where the upper (lower) sign is for B\\x (B\\y), respec- 
tively. 

From these simple first-order results, it is possible to 
understand the qualitatively different behavior obtained 
with local, single-ion, anisotropy from that obtained with 
global anisotropy. In the isotropic case Jj = OVj, the first 



FIG. 6: Plot at k B T/\J\ = 0.03 and J b /J = 0.1 of M/j versus 
r )B/\J\ for the AFM spin 1 dimer. Curves for B\\z (solid), 
B\\x (dashed), and the isotropic case (J b — 0, dotted) are 
shown. 



and second level crossings occur at — J and — 2 J, respec- 
tively. For each induction direction, the signs of the Jb 
and Jd contributions to B^^ + 2 J and B^^ + J are the 
same, whereas the signs of the J a and J e contributions 



to B. 



lc(l) 



2 J and B^{ 1] 



J are the opposite. 



In Figs. 6-15, we plot M/j and Cy/ks versus jB/\J\ 
for five low-T cases of AFM s\ = 1 dimers, taking 
UbT/\J\ = 0.03. The M/7 curves all exhibit the univer- 
sal features predicted by Eqs. (39) and (40). The corre- 
sponding Cy/kB curves also obey the universal features 
predicted in Eqs. (36)- (38). In these figures, only one of 
the five anisotropy interactions Jj is non-vanishing, and 
we take Jj/J = 0.1, for j = b, d, a, e, and c, respectively. 
Unlike the case of si = 1/2 dimers, for which J c merely 
renormalizes Jd, for si — 1 all of these interactions lead 
to distinct anisotropy effects in the low-T magnetization 
and specific heat. Of course, since the case J c = 0.1J 
(and the remaining Jj = 0) can be evaluated by setting 
Jd = 0.05 J and J e = — 0.05J for any si, it is not really 
distinct from the other anisotropy interactions, but only 
represents that particular combination of global and local 
azimuthal anisotropy interactions. 

We first show the results for the global anisotropy in- 
teractions. In Figs. 6 and 7, we plot M(B) and Cv(B) 
for J b = 0.1J (and all other Jj = 0) at k B T/\J\ = 0.03 
for the AFM si = 1 dimer. The level-crossing induc- 
tion increases monotonically with level-crossing number, 
with the largest induction required occurring for B\\z, 
as for si = 1/2, quantitatively consistent with Eqs. (54)- 
(57). In Figs. 8 and 9, the corresponding results for 
Jd = 0.1J (and the other Jj = 0) are shown. Again, 
the level-crossing induction increases monotonically with 
level-crossing number, and the largest induction required 
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k T/| J 1=0.03 




FIG. 7: Plot at k B T/\J\ = 0.03 and J 6 /J = 0.1 of C v /k B 
versus -yB/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (J b — 0, dotted) 
are shown. 




g^B/|J| 



FIG. 8: Plot at k B T/\J\ = 0.03 and J d /J = 0.1 of M/ 7 
versus ^/B/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (Jd = 0, dotted) 
are shown. 

for each level crossing occurs for B \ \ x , quantitatively con- 
sistent with Eqs. (54)-(57), as for si = 1/2. We note that 
the solid curves for B\\z are nearly identical with the 
dotted curves for the isotropic case Jj = OVj. However, 
there is a very small distinction between these M and 
Cy curves. This distinction is also present for s\ = 1/2, 
as shown in Figs. 3 and 5, and is a detail of the shapes 
of the M and Cy curves arising from a slightly higher 
effective temperature with Jd ^ than in the other cases 
studied. 

Next, in Figs. 10-13, we show the corresponding curves 



kJ/|J|=0.03 B ||x 




g^B/ui 



FIG. 9: Plot at k B T/\J\ = 0.03 and J d /J = 0.1 of C v /k B 
versus ^B/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (Jb = 0, dotted) 
are shown. 



for M(B) and Cy(B) at low T for the local anisotropics 
J a = 0.1J and J e = 0.1J, respectively, with the other 
Jj = 0. In contrast to the global anisotropics, the in- 
crease in B required with the level crossing is not mono- 
tonic in the level-crossing number, quantitatively consis- 
tent with Eqs. (54)- (57). In particular, we note that in 
each of these figures, the deviations in the level crossings 
B^i + 2 J and B^ + J are opposite in sign. 

A careful examination of the exact formula for Cy re- 
veals that it vanishes at precisely the level crossing induc- 
tions as T — > 0. Similarly, the magnetization curve at- 
tains its midpoint values \ and | between the two neigh- 
boring steps at the level crossing inductions as T — > 0. 

Finally, in Figs. 14 and 15, we show the corresponding 
M(B) and Cy(B) curves for J c = 0.1J, corresponding 
to the parameter choices J e = —0.05 J and Jd — 0.05J. 
With this combination, the magnetization and specific 
heat curves shown in Figs. 14 and 15 resemble those 
for J d = 0.1J shown in Figs. 8 and 9, except that the 
second magnetization steps and the second set of specific 
heat peaks have nearly the same behavior as do the first 
magnetization steps and specific heat double peaks. 

In short, the case si = 1, for which the exact ex- 
pressions from which the eigenvalues can be readily ob- 
tained are sufficiently short as to be writable on paper, 
are sufficient to exhibit the very different behaviors ob- 
tained from the single-ion, local spin anisotropy interac- 
tions from those obtained from the global spin anisotropy 
interactions. As s\ increases beyond 1, the situation be- 
comes not only more complicated, but also more inter- 
esting, as shown in the following. 
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FIG. 10: Plot at k B T/\J\ = 0.03 and J a /J = 0.1 of M/j 
versus ^B/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (J a = 0, dotted) 
are shown. 
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FIG. 12: Plot at k B T/\J\ = 0.03 and J e /J = 0.1 of M/j 
versus ^B/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (J e = 0, dotted) 
are shown. 




g^ R B/|J| 



FIG. 11: Plot at k B T/\J\ = 0.03 and J a /J = 0.1 of C v /k B 
versus ^B/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (Jb = 0, dotted) 
are shown. 



FIG. 13: Plot at k B T/\J\ = 0.03 and J b /J = 0.1 of C v /k B 
versus ^B/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (Jb = 0, dotted) 
are shown. 



EXACT NUMERICAL RESULTS FOR SPIN 5/2 

For si = 5/2, one of the cases of greatest experimen- 
tal interest, when Ti. a and H. e are present, none of the 
allowed s, m values is a true quantum number. That is, 
H a and H e cause all of the states with nominally odd or 
even s to mix with one another. For B\\i for i = x,y, z, 
this simplifies in the crystal representation as for si = 1, 
since only states with odd or even m in the appropriately 
chosen representation can mix. By using symbolic ma- 
nipulation software, it is possible to solve for the exact 



eigenvalues of the s\ — 5/2 dimer. However, because the 
analytic expressions for the eigenvalues are much more 
complicated than those for s\ = 1 presented in Appendix 
A, we shall not attempt to present them, but will instead 
focus upon their numerical evaluation for specific cases. 

To first order in the Jj, the first three level crossings 
for B\\i with i — x,y,z are 

7<L = -J+^-J b , (58) 

7<)U = -'-^-T*T ±1J T' (59) 
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FIG. 14: Plot at fe B T/|J| = 0.03 and J c /J = 0.1 of M/y 
versus yB/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (J c — 0, dotted) 
are shown. 
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FIG. 15: Plot at k B T/\J\ = 0.03 and J c /J = 0.1 of C' v /k B 
versus yB/\J\ for the AFM spin 1 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (Jb = 0, dotted) 
are shown. 
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In Figs. 16-23, we plot M/7 and Cy /k B versus jB/\J\ 
for four low-T cases of AFM s\ = 5/2 dimers, = 0.1 J, 
Jd = 0.1J, J a = 0.1J, and J c = 0.1J, respectively, and 



FIG. 16: Plot at k B T/\J\ = 0.03 of M/7 versus jB/\J\ for 
the AFM spin 5/2 dimer with J b /J = 0.1. B\\z (solid), 
(dashed), and the isotropic case (Jb = 0, dotted) are shown. 
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FIG. 17: Plot at k B T/\J\ = 0.03 and J b /J = 0.1 of C v /k B 
versus jB/\J\ for the AFM spin 5/2 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (J b — 0, dotted) 
are shown. 



the other Jj = 0, taking k B T/\J\ = 0.03. In Fig. 24, ex- 
amples of M(6) at fixed B and <j> = are shown. Figures 
16-24 are sufficient to distinguish the more interesting 
local spin anisotropy effects in AFM dimers with higher 
si values from the non-existent or less interesting ones 
present with si = 1/2, 1, respectively. In Figs. 16, 18, 20, 
and 22, the solid and dashed curves represent the cases 
of B\\z and B\\x, and the dotted curve is the isotropic 
case, Jj = OVj, as in Figs. 2-15. Because of the num- 
ber of peaks in the specific heat curves, in Figs. 17, 19, 
21, and 23, we ol=nly showed Cy jk B versus jB/\J\ for 
B\\z (solid) and B\\x (dashed). The isotropic curve with 
Jj = L\/j was published previously. [14] That curve has 
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FIG. 18: Plot of M/7 versus -)B/\J\ for the AFM spin 5/2 
dimer at k B T/\J\ = 0.03 with J d /J = 0.1. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (Jb = 0, dotted) 
are shown. 
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FIG. 19: Plot at k B T/\J\ = 0.03 and J d /J = 0.1 of C v /k B 
versus -yB/\J\ at k B T/\J\ = 0.03 for the AFM spin 5/2 dimer. 
Curves for B\\z (solid), B\\x (dashed), and the isotropic case 
(Jb — 0, dotted) are shown. 



minima at jB/\J\ = s for s = 1, . . . , 5, each of which is 
central to double peaks. H a and H b are invariant under 
x <-> y, so curves for in Figs. 16, 17, 20, and 21 

are identical to those for Since Ttd and H c are odd 

under x <-> y, the curves in Figs. 18, 19, 22, and 23 
correspond to 23||y with J^/J = — 0.1J, J c /J = —0.1, 
respectively. 

We first examine the global spin anisotropy effects of 
Hb and Hd in Figs. 16-19. These figures exhibit the 
same behavior shown for si = 1/2, 1 in the correspond- 
ing Figs. 2-9. Note that 5' c 5/2 (6>, 0) is largest for 6 = tt/2 
with J b = 0.1J and for 9 = with Jd = 0.1J, and in- 



creases monotonically with s, as for s\ = 1 in both cases, 
nearly quantitatively consistent with Eqs. (58)- (63). By 
contrast, B l s c 52 (6,0) for 9 = tt/2 with J b = 0.1 J and for 
8 = for Jd = 0.1J are nearly indistinguishable from the 
isotropic case, also nearly quantitatively consistent with 
Eqs. (58)-(63). 

In contrast, the local field anisotropy interactions show 
very different and much more interesting behaviors. In 
Figs. 20-21, we present our results for the effects of the 
single-ion axial anisotropy interaction 7i a , Eq. (5). As 
in Fig. 10 for si = 1, B| c y 2 (fl, 0) + sJ changes sign 
with increasing s. -E>i c 5 / 2 (0, 0) + J < 0, whereas for 
s > 2, B^ 5 , 2 (0, 0) + sJ > and increases monotoni- 
cally with s. For 9 = ir/2, nearly the opposite situa- 
tion occurs. i?' c 5 y 2 (7r/2, 0) + sJ is positive for s = 1 
and decreases monotonically with increasing s. In both 
cases, \B l s 5 / 2 (@i 0) + i s a minimum for s = 2. These 
local axial anisotropy effects, consistent with Eqs. (58)- 
(63), are very different than the global anisotropy ones 
pictured in Figs. 16-19. They are also much richer and 
interesting than the corresponding case for s\ = 1 pic- 
tured in Figs. 10 and 11. 

In Figs. 22-23, we present our M/7 versus jB/\J\ 
results for the case of the local azimuthally anisotropic 
exchange interaction, H c , Eq. (8), evaluated for AFM 
dimers at ksT/\J\ = 0.03 with J c = 0.1J and the re- 



maining Jj = 0. -835/2 
Eqs. (58)-(63) by setting J d , J, 



for this case is obtained from 
±J c /2, respectively, 
and J a = J b = 0. -B' c 5 / 2 (0,0) is nearly indistinguishable 
from the isotropic case, as in Figs. 12 and 13 for si = 1, 
except for some minor curve shape effects far from the 
step midpoints. B\ c 5 ^ 2 (tt/2 7 0) + sJ is always positive, as 
for Si = 1 shown in Figs. 12 and 13, but has a minimum 
at s = 3. This is also in stark contrast to the monotonic 
global anisotropy behavior for si = 5/2 seen in Figs. 
16-19. Both of these behaviors are nearly quantitatively 
consistent with Eqs. (58)-(63) as modified to include J c . 

In contrast, the local field anisotropy interactions show 
very different behaviors. In Figs. 20 and 21, our results 
for the effects of the single-ion axial anisotropy interac- 
tion H a , Eq. (5), are shown. As in Figs. 10 and 11, 
the induction anisotropy effects of 7i a change sign with 
increasing B. In Fig. 20, for B\ \z the first magnetization 
step appears at a lower value of \B\ than in the isotropic 
case, and for steps 3-5, there is a monotonic increase in 
the extra field required for each step. The opposite is true 
for B\\x, for which the first step appears at a larger |B| 
value than for the isotropic case, and subsequent magne- 
tization steps appear at monotonically decreasing values 
of \B\. A cross-over occurs at about the second step, 
for which the effects of this type of anisotropy are small. 
The corresponding shifts in the positions of the double 
peaks in the specific heat are shown in Fig. 21. These 
local axial anisotropy effects are very different than the 
global anisotropy ones pictured in Figs. 16-19. 
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FIG. 20: Plot at k B T/\J\ = 0.03 and J a /J = 0.1 of M/ 7 
versus 7 B/|J| for the AFM spin 5/2 dimer. Curves for B\\z 
(solid), B\\x (dashed), and the isotropic case (J a = 0, dotted) 
are shown. 
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FIG. 22: Plot at fc B T/|J| = 0.03 of M/ 7 versus 7 B/| J| for 
the AFM spin 5/2 dimer with J c /J = 0.1. B\\z (solid), Bp 
(dashed), and the isotropic case (Jb = 0, dotted) are shown. 
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FIG. 21: Plot at k B T/\J\ = 0.03 and J a /J = 0.1 of Cv/k B 
versus 7 B/| J\ at k B T/\J\ = 0.03 for the AFM spin 5/2 dimer. 
Curves for B\\z (solid), B\\x (dashed), and the isotropic case 
(Jb — 0, dotted) are shown. 



FIG. 23: Plot at k B T/\J\ = 0.03 and J c /J = 0.1 of C v /k B 
versus -yB/\J\ at k B T/\J\ = 0.03 for the AFM spin 5/2 dimer. 
Curves for B\\z (solid), B\\x (dashed), and the isotropic case 
(Jb — 0, dotted) are shown. 



In Figs. 22 and 23, we present our M and Cy results 
for the case of the local azimuthally anisotropic exchange 
interaction, H c , Eq. (8), evaluated for AFM dimers at 
UbT/\J\ = 0.03 with J c = 0.1J, or equivalently with 
Jd = 0.05 J and J e = — 0.05 J (and the remaining Jj = 0). 
For B\\z, there is almost no change from the isotropic 
occurred with the global azimuthal anisotropy 
pictured in Figs. 128 and 19. For the field required 

for each step is larger than in the isotropic interaction 
case, as in Figs. 14 and 15 for si = 1, and for global az- 
imuthal anisotropy shown in Figs. 18 and 19. However, 
in this case, the extra induction required for each level 



crossing is non-monotonic in the crossing number, with 
the largest extra induction required for the first crossing, 
and the minimum extra induction required for the inter- 
mediate, third crossing. This is in stark contrast to the 
monotonic global anisotropy behavior for Si = 5/2 seen 
in Figs. 16-19. Although not pictured explicitly, the be- 
havior for J e = 0.1J with the other Jj — is rather like 
that of the J a — 0.1J curves pictured in Figs. 20 and 21 
with x <-> z, differing in ways similar to those differences 
between the J a = 0.1J and J e = 0.1J curves pictured for 
si = 1 in Figs. 10-13. 

Although not pictured explicitly, the behavior for J e = 
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FIG. 24: Plot of M/7 at <^ = versus 9/n near the second step 
at jB/\J\ = 2.0,2.1 and k B T/\J\ = 0.03 for each Jj/J = 0.1 
with i — a,b, c, d. Curves are labelled with j, jB/\ J\ values. 

0.1J with the other Jj = is rather like that of the 
J a = 0.1J curves pictured in Figs. 20 and 21 with z <-> 
x, differing in ways similar to those differences between 
the J a — 0.1J and J e = 0.1J curves pictured for si = 
1 in Figs. 10-13. As indicated in Eqs. (58)-(63), the 
£^ c 5 / 2 (0, 0) are nearly independent of J e . However, the 
B' c 5 ^ 2 (7r/2, 0) + sJ with J e = 0.1J are nearly three times 
as large as for J a = 0.1J case, so that the minimum 
\B 1 ^ 5 ^ 2 (n/2, 0) + sJ\ is also a minimum for s — 2. More 
details are given in Section VII and Appendix D. 

In addition, the angular dependencies of M are differ- 
ent for each of the four Jj we presented for si — 5/2. In 
Fig. 24, we present the results for |M(£?,6>, 4> = 0)\/"f 
versus 9/tt near the second level crossing at -f\B/J\ = 
2.0, 2.1 and k B T/\J\ = 0.03 for each of the four si = 5/2 
AFM magnetization cases pictured in Figs. 16, 18, 20, 
and 22. Note that \M\(B, tt - 9, 0)| = \M(B,8,0)\. 
The J a and J curves, while rather similar at j\B/J\ = 
2.0, are very different at -f\B/J\ = 2.1. Hence, 
M(B) depends strongly upon the particular type of spin 
anisotropy. 



ANALYTIC RESULTS FOR WEAKLY 
ANISOTROPIC DIMERS OF ARBITRARY SPIN 

Induction representation eigenstates first order in 
the anisotropies 

Since the diagonalization of the Hamiltonian matrix 
is difficult for an arbitrary magnetic field H direction 
and for an arbitrary combination of spin anisotropy in- 
teractions, and must be done for each value of s\ sepa- 
rately, it is useful to consider a perturbative solution in 
the relative strengths Jj/J of the the anisotropy interac- 
tions. We nominally assume \ Jj/J\ <C 1 for j = a,b,d,e. 
However, to compare with low-T M(B) and Cv(B) ex- 



periments at various H directions and magnitudes, one 
cannot take B to be small. In order to incorporate an ar- 
bitrary B accurately, we rotate the crystal axes (x, y, z) 
to (x' ,y' , z'), so that B — Bz' . The rotation matrix 
and a brief discussion of its ramifications are given in 
Appendix B. 

In these rotated coordinates, the Zeeman interaction 
~jBS z > is diagonal. We therefore denote this represen- 
tation as the induction representation. The Hamiltonian 
Tt' in this representation is given in Appendix B. In the 
induction representation, we choose the quantum states 
to be |</?™). In the absence of the four anisotropy inter- 
actions Jj, Tt' = Tt' is diagonal, 

WolvD = K' m \v7), (64) 

where 

£™' (0) = -Js(s + l)/2-jBm. (65) 

The operations of the remaining terms in Tt' on the 
eigenstates Yp™) are given in Appendix C. The first order 
correction to the energy in this representation is E™'^ — 
(ip™\H'\ip™), which is found to be 

= -y[2*(* + l)-l]-Y[*(« + l)-l] 

+ ^-[m 2 + s(s + 1)-1] 
+ i[ S (.s + l)-3m 2 ] 

x (jbM cos2 9 + Jd^ sin2 6 cos(20)) , (66) 

where 

•Kfa = Jb + Us, Sl Ja, (67) 
= Jd + Ct s , Sl J e , (68) 

and a s , Sl is given by Eq. (29). 

Since the 0, </> dependence of ET arises from the 
term proportional to J^'* 1 cos 2 9 + sin 2 9 cos(20), 
it is tempting to think that the thermodynamics with 
Jde 1 = and B\\z are equivalent to those with J^ 1 = 
and B\\x. However, as shown explicitly in the following, 
the 9, ^-independent parts of Eq. (66) strongly break this 
apparent equivalence, causing the Bl c si (9,(j)) for those 
two cases to differ. This implies that J and Jd are in- 
equivalent, as are J a and J e , even to first order in the 
anisotropy strengths. 

First order thermodynamics 

To first order in the anisotropy interactions, s and m 
are still good quantum numbers, so the partition function 

^EE e-'^ (69) 

s— m— — s 
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where E" 



E' s 



m,(0) 



p m,(l) 



Although it is difficult to 



perform the summation over the m values analytically, it 
is nevertheless elementary to evaluate Z numerically for 
an arbitrary B, 9, (ft, and T from the eigenstate energies. 
The magnetization is obtained from 



M(B,6,< 



2si 
s=0 m= — i 



7 V 



me 



(70) 



Similarly, the specific heat to first order in the 
anisotropy interactions is found from 



2s\ S 



s—0 m— - 

2si s 



"(£ E K,s,e-^) }. (71) 

s— m— — s 

As a test of the accuracy of this first-order calculation, 
we have compared the first-order and exact M(B) ob- 
tained for the si = 5/2 dimer with Jd = 0.1J and B\\z 
in Fig. 25. The corresponding comparison between the 
first-order and exact Cy(B) is shown in Fig. 26. We 
see that the curves evaluated using the first-order and 
the exact expressions for M and Cy with si — 5/2 are 
indistinguishable at fceT/l J| = 0.03. The Cy curves are 
noticeably different at k B T/\J\ = 0.1 for jB/\J\ < 0.4, 
and at UbT/\J\ — 0.3, they are noticeably different for 
jB/\J\ < 2.6. Corresponding noticeable differences in 
the M curves at the same B values appear at T values 
roughly three times as high as in the Cy curves. 

At very low T, fcsT/|J| <C 1, the most important 
states in this perturbative scheme are the minima for 
each s value, El Si , which determine the first-order level 
crossings. As T — > 0, we can ignore all of the m ^ s 
states in Eqs. (69)-(71). 



Level crossings first order in the anisotropics 
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FIG. 25: Comparison of M/7 versus yB/\J\ obtained using 
the first-order asymptotic form (dotted) with the exact cal- 
culation (solid), for the si = 5/2 AFM dimer with Jd = 0.1J, 
J a = J b = Je = 0, at k B T/\J\ = 0.03,0.1,0.3,0.8, as indi- 
cated. Inset: expanded view of the region < 7-B/l J\ < 1.5. 
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We can find an expression for the s th AFM level cross- 
ing at the induction B l s c ^ to first order in the anisotropy 
interactions for a general si spin dimer by equating 

Et^+EtiiUoE:zi' {o) +E s -iii\ 

IB 1 *'? = -Js-J b /2-c SiSl J a 

- ^ S ~ 3 ^ [Jb cos 2 9 + J d sin 2 9 cos(20)] 

+3c s , Sl [J a cos 2 9 + J e sin 2 9 cos(2</>)], (72) 
[3 + 3s - 5s 2 - 4s 3 + 4si(si + 1)] 



2(2s+ l)(2s + 3) 



(73) 



This expression is consistent with those obtained for 
si = 1/2,1 given by Eqs. (52) and (54)-(57). In ad- 
dition, this expression is nearly quantitatively in agree- 
ment with the M(B) and Cy(B) behaviors pictured for 



FIG. 26: Comparison of Cy/fcs versus 7B/I J| obtained using 
the first-order asymptotic form (dotted) with the exact cal- 
culation (solid), for the si = 5/2 AFM dimer with J d = 0.1J, 
Ja = J b = Je= 0, at k B T/\J\ = 0.03,0.1,0.3, as indicated. 



si = 5/2 in Figs. 16-23. We note that jB 1 ^ 1 ^ contains 
the 9, (^-independent terms, — Js — Jt/2 — c s , Sl J a , which 
distinguish Jf, from Jd and J a from J e . 

In particular, we note that the single-ion anisotropy 
interactions behave very differently with increasing step 
number than do the global anisotropy interactions, espe- 
cially for large si. For the three cases we studied in de- 
tail, for si = 1/2, c : 1/2 = 0, so that the local anisotropy 
terms are irrelevant, for si = 1, ci t i — g and C2.1 = — \ 
have different signs, and for s\ — 5/2 as in Fe2 dimers, the 
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first three c s ^/2 coefficients are y|, — and — 1|, respec- 
tively, the second being an order of magnitude smaller 
than the other two, and opposite in sign from the first. 
For si — 9/2 dimers such as [Mn^, the first four c s .g/2 
are Tp, |, — |, and — §|, which changes sign between 
s = 2 and s = 3, where its magnitude is a minimum. 

This is in sharp contrast to the global anisotropy inter- 
actions, for which the analogous coefficient (4s — 3)/2 in- 
creases monotonically from 5 to | as s increases from 1 to 
2, independent of s\. These differences should be possi- 
ble to verify experimentally in careful low-T experiments 
at high magnetic fields applied at various directions on 
single crystals of those s\ — 5/2 Fe2 and si = 9/2 [Mn^ 
dimers for which | J| is sufficiently small. 



Level crossings to second order in the anisotropy 
energies 

To aid in the analysis of experimental data, we have 
extended this perturbative calculation to second order. 
Since we expect the \Jj/J\ < 0.1 in most circumstances, 
this extension ought to be sufficient to accurately analyze 
most experimentally important samples. To second or- 
der in the anisotropy interactions, the cigenstate energies 
E™g[ 2 ^ are given in Appendix C. We note that the i?™i 1 2 ' ) 
contains divergences at jB/\J\ = 0, 2s— 1, 2s + 3, s— 1/2, 
and s + 3/2, so that near to those values, one would need 
to modify the perturbation expansion to take proper ac- 
count of the degeneracies. Hence, the expressions for 
E™sn cannot be used in the asymptotic expressions for 
the thermodynamics, Eqs. (69)-(71). However, as the sth 
AFM level crossing occurs approximately at 7-B/l J\ = s, 
which is far from any divergences, we can safely use this 
second order expansion to obtain an expression for the 
level crossings second order in the anisotropy interaction 
energies. We find 



-Ja/i 



(74) 



where jB^i^ is given by Eq. (73). 

The full expression for Efcg is given in Appendix D. 
From this expression, it is easy to see that for si = 1/2, 

^11/2 = 0) + \fi{®i 0)> wri crc fi and / 4 are given 

in Appendix D. For B\\z, E^^L — -ijJ/|J|, and with 

B\\x,y, -E^ ^ = ^ J<l) 2 l\J\i in agreement with 

the expansion to second order in the Jj of our exact for- 
mulas in Eq. (52). However, the second order functions 
have more complicated 0, <j) dependencies than do the first 
order B^s} in Eq. (73). We have also explicitly checked 
each formula in Appendix D for si = 1 and s = 1,2. 
Thus, Eq. (74) is a highly accurate expression for the 



full 9, 4> dependencies of all s = 1, . . . , 2s\ + 1 level cross- 
ings of a single crystal dimcr of single ion spin s\. 

By superposing this non-universal level-crossing for- 
mula, Eq. (74), combined with the universal behavior 
presented in Eqs. (??)-(40), it is easy to use our results 
in accurate fits to experimental data at low temperatures 
and high magnetic fields on dimers with arbitrary si val- 
ues. 



SUMMARY AND CONCLUSIONS 

In summary, we solved for the low-temperature mag- 
netization and specific heat of equal spin si antiferro- 
magnetic dimcr single molecule magnets, including the 
most general forms of anisotropic spin exchange inter- 
actions quadratic in the spin operators. The magneti- 
zation and specific heat exhibit steps and zeroes, respec- 
tively, at the non-universal level-crossing induction values 
B 1 g Si (0, </>), but the magnetization steps and their mid- 
point slopes, plus the two peaks surrounding the specific 
heat zeroes all exhibit universal behavior at sufficiently 
low temperatures. Local (or single-ion) anisotropy in- 
teractions lead to low-temperature magnetization step 
plateaus that have a much richer variation with the mag- 
netic induction B than do those obtained from global 
anisotropy interactions, provided that si > 1/2. We de- 
rived simple, accurate asymptotic analytic expressions 
for the low-temperature magnetization and specific heat 
for the most general quadratic anisotropic spin interac- 
tions at an arbitrary B, and an accurate expression for 
Bg C sl (6, </>), enabling fast and accurate fits to experimen- 
tal data. 

There were two low-T M(B) studies of Fe2 dimers. [16, 
17] For / u-oxalatotetrakis(acetylacetonato)Fe2, all five 
peaks in dM/dH were measured in pulsed mag- 
netic fields H. These evenly spaced peaks indi- 
cated little, if any, spin anisotropy effects. [16] On the 
other hand, studies of the first 2-3 dM/dH peaks 
in powdered samples of [Fe(salen)Cl]2, where salen is 
N, A r '-ethylenebis(salicylidcnciminato), were much more 
interesting. [17] These data showed a broad first peak at 
B = 17 — 20T that was only partially resolvable into 
two separate peaks, followed by a sharp second peak 
at B = 36 T, consistent with local axial anisotropy of 
strength \J a /J\ w 0.1, as obtained from the derivatives 
of the curves shown in Fig. 20 From Eqs. (58)-(63), 
one could also have |(J a ±3J3)/J|«0.1. These values 
might perhaps be combined with a smaller |J C /J|, ob- 
tained from the derivatives of the curves pictured in Fig. 
22. 

Without a detailed single crystal study with the mag- 
netic field directed along a number of different crys- 
tal directions, it is impossible to determine the relative 
amounts of J a and J e that would best fit the data. How- 
ever, the existing data on [Fe(salen)Cl]2 appear to be 
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inconsistent with a predominant global anisotropy inter- 
action of either type, as obtained from the derivatives of 
the curves shown in Figs. 16 and 18. Only single crystal 
studies could determine if a small amount of such global 
anisotropy interactions were present in addition to one 
or more presumably larger local spin anisotropy interac- 
tions. In comparing the two materials cited above, it ap- 
pears that the interaction of a Cl~ ion neighboring each 
Fe 3+ ion leads to strong local (or single-ion) anisotropy 
effects. In order to verify this hypothesis and to eluci- 
date the details of the interactions, further experiments 
using single crystals in different field orientations on this 
and related Fe2 dimers with 1-3 similarly bonded Cl~ 
ions are urged. [19, 20] We also urge single crystal data 
on some of the si — 9/2 [Mn^ dimers, [21, 22] as well 
as on si = 1/2 dimers lacking in predicted local spin 
anisotropy effects. [11-14] To aid in the fits, we derived 
simple, useful formulas for the magnetization and specific 
heat at low temperature and sufficiently large magnetic 
induction, and accurate formulas for the level crossing 
inductions. 

With local anisotropy interactions, the total spin s is 
not a good quantum number, potentially modifying our 
understanding of quantum tunneling processes. It might 
also be possible to fit a variety of experimental results us- 
ing a smaller, consistent set of model parameters. [6] We 
emphasize that the study of SMM dimers, for which the 
most general anisotropic quadratic exchange interactions 
can be solved exactly, may be our best hope for attain- 
ing a more fundamental understanding of the underlying 
physics of single molecule magnets. 
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APPENDIX A 
Specific heat details for si = 1/2 

We first present the numerators of the exact expres- 
sions for the specific heat with s\ = 1/2 and B\\i for 
i = x,y, z. We have 

K, v = Fly + \{J + 2J y , x ) 2 eW J ^- J *^ 

+F x , y smh((3F x , y ) e - 0J *-v 

x [(J + J x ,y)e~P J + (J x , y - 2J v , x )e 20J y'*] 



+ ^e- 0J *-y cosh((3F x , v ) 

x([(J + J^) 2 + ^ 2 Je-^ 

+ [{2J y , x -J x , v Y+Fiy^*), (75) 

N z = F z 2 cosh(2/?F z ) + ^ e -' 3 ( J + 2J ») 
+ icosh(/3iy (a z (J 2 +F 2 ) 
+e -0(J+J»)j(j + 2J 6 )) 
+F Z sinh(/3F z ) ( J b A z + J e -^ J+J ^ . (76) 

Eigenvalues for si = 1 

In the remainder of this appendix, we provide the de- 
tails of our exact results for si = 1. The cubic equation 
for the three s = 1 eigenvalues is given by 

e„ = -J-Jt-Ja + K, for n = 2, 3,4, (77) 

= -\ 3 n -(Ja- J b )\ 2 n +K[b 2 + (J d - Jef] 

+(J a - J b )[b 2 z + (J d - J e f] 

-(J d -J e )(b 2 x -b 2 y ), (78) 

where the bi = jBi as in Eq. (??). For B\\z, the cubic 
equation is easily solved to yield 

A„ = J„- J a ,±V^ + (J d - J e )2. (79) 

For B\\x, y, we have 

A„ = ±{J e -J d ),-J y , x ± ^V+J 2 XtV , (80) 

Jx,y — -j[Ja — JbT {Jd - Je)], (81) 

where J x (J y ) corresponds to the upper (lower) sign. The 
rank ordering of these eigenvalues depends upon B and 
the various anisotropy parameters, of course. 

The six eigenvalues for the mixed s — 0, 2 states satisfy 

e n - -3J-^J a -2J h + X n (82) 
We then define 

J a = —^-Ja, (83) 

Jb - Jb + Ja/3, (84) 
Jd = Jd + Je/3, (85) 

Je = -^Je, (86) 
J = J- \ja, (87) 
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and obtain the Hcrmitian matrix M for the six s = 0, 2 
states, the eigenvalues of which are the A„. For brevity, 
we let 



np 

"in 


T7 Tr 1 Tl/l ('• (~\ G /J 

IVU\y \ JJU LUD (/ , 






— L> bill (76 , 




b 3 




(90) 


a 




(91) 


d 3 


= 3Jd, 


(92) 


d 6 


= \Z6Jd, 


(93) 


e 


= Je- 


(94) 



CO 



-2 J 2 (J 2 + 9J 2 ) + 8J 2 (2J 2 + 6J 2 + J 2 ) 

+&V&j a j b j d j^j 

-2J d (6(6 2 -3J 2 + J 2 )(6 2 -6 2 ) 
-V6 JaMb 2 - 3b 2 z - 18 Jj)] 
-12 J ( J b b 2 (b 2 -3b 2 z ) 



(102) 



-ZJd{b 2 
-Jl{b 2 ^ 



-Jl 
3b 2 z - 

2 



2 J b 2 + 12J, 2 



Then, the matrix M is given by 



M = 





-b± 


-d 6 








-e \ 


-61 


Qi 1 


-b 3 


-d 3 








-de 


-b% 


Ql 


-b 3 


-d 6 


—a 





-d 3 


-b% 


Q\ 


-b± 











-d e 


-bl 


Q 2 -2 


— e 


V ~e 





—a 





— e 


Q° + 3 J / 



.(95) 



After some rearrangement, the A„ are found to satisfy 

6 



-3J b J^(7b z - 3fo 2 - 18 J d 
-J 2 a ((b 2 -3b 2 z + 2J 2 ) 2 + 36J 2 (b 2 z -J 2 )) 
-4J 2 (2fe 2 (fe 2 - 36 2 + J 2 ) - 2J 2 (2 J 2 + J 2 - b 2 - 3b 2 z 
+6J 2 d (-7b 2 + 36 2 + 4J 2 + 18J 2 + 3J 2 ) 

+ (6J b J d - VGJaJeM -6 2 )) 

+24jJ d (6 2 -6 2 )(6 2 -3J 2 -2J e 2 ) 
+2V6J a J e (b 2 x - b 2 y )(b 2 - 36 2 - 6J 2 ) 
-6J 2 (& 2 -6 2 ) 2 

+8V6J a J b J d J e (2b 2 - 3b 2 z + J 2 - 9Jj). 



(103) 



= ^c P (A„f, 

C6 = 1, 

c 5 = -2J b -3J, 

c 4 = - J 2 - 7J 2 - 21 J 2 - 2J e 2 - 56 2 , 

c 3 = 3J(7J b 2 + 21 J, 2 + 5b 2 ) 

+2 J b (-J 2 + 8 J 2 + 12 J 2 + 2 J 2 ) 
+4V6J a J d J e + 21 J d (fe 2 - 6 2 ) 
+3J b (6 2 + 76 2 ), 



(96) 

(97) 
(98) 
(99) 



B along a crystal axis 



We recall that the six Hamiltonian matrix eigenvalues 
e n are generally given by Eq. (82). For the special case 
B\\z, M is block diagonal. The resulting eigenvalues of 
the nominal s — 2, m — ±1 notation are 



A 



(100) 



6.8 



- J b + J6 2 + 9J 2 



(104) 



C2 



3J~(j b [- 



2 J 2 + 18J 2 + 7(fe 2 - 3& 2 )] 



The remaining four states are obtained from 



Cl 



-21 J> 2 - b 2 y ) 

+J 2 (3J 2 + 9J 2 + 2fo 2 + 36 2 ) 
+2J 2 [4J 2 + 3J 2 + ll(6 2 - 36 2 )] 
+2J 2 [54(J 2 + J 2 )+9J e 2 -12(6 2 - 
+2J e 2 (46 2 -36 2 ) + 46 4 
+2(3 J 6 J d - V6J a Je)(b 2 x -b 2 y ), 

-12 J~(V + 3J 4 + 18J 2 J 2 + 27 J 4 

3J 2 )(6 2 -36 2 ) 



36 2 )] 



(101) 



+2(J 6 2 



+ 12J b J d (b 2 x - b 2 ) 
-2J b (2(9J 2 + b 2 ){b 2 






4 


(105) 




= 1, 


(106) 




= -3J, 


(107) 


A; 2 


= -6 JJ b - J 2 - 8 J 2 - 12 J 2 - 2 J 2 


-4fe 2 ,(108) 


*? 


- 12J> 2 + J 2 +3J 2 ) 






+4J b (4fe 2 - J 2 )+4\/6J a J d J e , 


(109) 




- 24JJ b (-fo 2 + J 2 + 3J 2 ) + 4J 2 (6 2 






+8J 2 (— 26 2 + 2J 2 + 6J d + J 2 ) 






+8V6J a J&JdJ e - 


(110) 



Zb 2 z ) 



-J 2 (6 2 -36 2 )+6J 2 (fe 2 -56 2 ) 
+ (42J b J d + V6J Q J e )(6 2 -6 2 ) 



Similarly for B\\x,y, M is also block diagonal. The 
resulting eigenvalues of the nominal s = 2, m = ±1 no- 
tation are 



3 - 

\x,y _ j 

A 6,8 — 2 



2J b + + J^, 



(111) 
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Jx,y = -j{Jb±Jd)- 

The remaining four states are obtained from 



(112) 



n 



= 



x,y 



x,y 



x,y 

n ,y = 



x.y 

n-, y = 



(113) 

(114) 
(115) 

(116) 



x,y 



j2^ P - v (K' v r, 

p=0 

1, 

-3(J + J b T Jd), 
- J 2 - 2 J 2 - 12 J 2 - 2 J 2 + 3JJ b 
-4b 2 T 3J d (2 J b + 3 J), 
6J e 2 (J 6 T Jd) + 4fe 2 (3J + J 6 ± 3J d ) 
+ 12(J 2 + 3J 2 )(J + J b T Ji) 
- J 2 ( J & ± 3 J d )+ 4^6 J Q J d J e , 
b 2 [J 2 + 12 J(J b t 3J d ) + 8 J 2 + 6 J 2 
T24 J 6 J d T 2V6 J a J e ] 
+2J 2 [J a 2 - 4 J 2 - 6 J Jfc - 12 J 2 - 2 J 2 ] 
±12J d (2J b + 3J)(J 2 +3J d 2 ) 
T6J b Jd(Ja±6JJd-2J 2 ) 
-4V6J a JdJe(Jb T 3J d ). (118) 

Simple special cases 



When only one of the Jj 7^ 0, the eigenvalues for B\\z 
simplify considerably. For J b ^ 0, we have 

\ 7 n = 3 J + 2 J fc , 2 J 6 , J 6 ± 6, -2 J h ± 26. (119) 

For J d 7^ 0, we find 



A* = 0,3J,±y& 2 + 9J 2 ,±2^/& 2 + 3J 2 . (120) 
For J a 7^ 0, the eigenvalues can also be found analytically, 



We note that the cubic equation must be solved to obtain 
the ground state energy, since the s = m — state mixes 
with the s = 2, ±2 states. For B\\x, y, with J b ^= 0,J d ^ 
0, the eigenvalues are 



K' v = 3J + 2J b ,-2J b + J X)V ±^b 2 + Jl y ,(126) 

where J x . y is given by Eq. (112), and the other three 
eigenvalues satisfy the cubic equation, 

= -x 3 + (J b T3J d )x 2 +A(b 2 + J b 2 + 3Ji)x 

+4(b 2 -J 2 -3J 2 )(J b T3Jd). (127) 



(117) For J a 7^ 0, two of the eigenvalues satisfy 



K' y = -y±^ 2 + m 



(128) 



and the other four eigenvalues satisfy the quartic equa- 
tion, 

= 8lx 4 -27(J a + 9J)x 3 -9(m 2 -9JJ a + l2J 2 )x 2 
+12(- Jl + 9JJ 2 - 9b 2 J a + 816 2 J)x 
+4(4J^ - 9JJ 3 + 18b 2 J 2 + 81b 2 JJ a ). (129) 

Finally, for J e ^ 0, two eigenvalues are 



Je 



K = - 2 ^±2b±±b, J -f + ^± s J\j> + Jl-jj a . 

(121) 

We note that the ground state energy in this case is 



El = ~J -Ja-\j\j 2 + Jl - JJa, (122) 

which explicitly involves mixing of the s = m = and the 
s = 2, m = states. For J e 7^ 0, three of the eigenvalues 
are 



Ki V = ±y + Vb 2 + J 2 /4,±y - Vb 2 + J 2 /4, 

(130) 

and the other four eigenvalues satisfy the quartic equa- 
tion, 

= x A - (3 J =F Je)x 3 - (46 2 + 4J 2 ± 3JJ e )x 2 
+4[3fc 2 J + J J 2 ± J e (b 2 - J 2 )]x 
±4J e (-3b 2 J±2b 2 J e + JJ 2 ). (131) 

At the first level crossing with B\\z, we have for J a ^ 



jB\ c f = i(j+[9J 2 + 4J 2 -4JJ a ] 1 /2). (132) 

At the second level crossing, simple formulas are only 
obtained for B\\z with one Jj 7^ 0. For B\\z and J a 7^ 0, 
J b 7^ and J^ 7^ 0, respectively, we have 



(123) 
(124) 



and the remaining three satisfy the cubic equation, 

= x 3 - 3Jx 2 - Ax(b 2 + J 2 ) + 4J( J 2 + 3fc 2 ).(125) 



-2J-J a , 
-2 J — 3 Jf,, 



(133) 

1/2 



i (20 J 2 - 30J 2 + 8[4J 4 - 6 J 2 J 2 ] 1 ^ . 
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APPENDIX B 



Global Hamiltonian matrix elements 



Rotation to the induction representation 

The rotation from the crystal representation to the in- 
duction representation is obtained from 



The operations of the rotated global anisotropy inter- 
actions upon these states may be written as 



ni\v>T) 



Jb 

' 4 



(4m 2 + 2[s{s + 1) - 3m 2 ] sin 2 6») \<p?) 



-sin(20) (2m + cr)AJ">™+ CT ) 



x 
V 



cos 


9 cos <fi 


— sin<^ 


sin 9 cos 4> \ 


<"\ 




COS 


9 sin <f> 


COS0 


sin sin <f> 




+ sin 2 




su\9 





cos 9 / 







<T=±1 



+ sin 2 9 J2 F s m W7 +2a ) 



<T=±1 



(137) 



(134) 



leading to B = Bz'. [25] This operation is equivalent to a 
rotation by — ir/2 about the z axis, a rotation by about 
the transformed x axis, and then a rotation by ir/2 — <f> 
about the transformed z axis. [26] In effect, in using the 
above rotation matrix, we made the arbitrary choice that 
the rotated z axis lies in the x'z' plane. After the above 
rotation, it is still possible to rotate the crystal by an ar- 
bitrary angle \ about the z' axis, keeping B\\z' . Hence, 
there are in effect an infinite number of equivalent rota- 
tions leading to B = Bz' . The resulting Hamiltonian 
matrix will then have off-diagonal elements that depend 
upon x- However, all such rotations necessarily lead to 
the identical, x-independent, set of eigenvalues of the re- 
sulting diagonalizcd Hamiltonian matrix. We have ex- 
plicitly checked that the above rotation gives the exact 
cubic expression, Eq. (43), for the s\ = 1/2 eigenvalues, 
and also leads to the correct eigenstate energies second 
order in each of the Jj for si = 1. We also showed explic- 
itly that x does not enter the eigenstate energies second 
order in J for arbitrary s,si,m. 



Global Hamiltonian 



The global axial and azimuthal anisotropy interactions 
in the rotated frame are 



Hi 



-J b [s 2 z , cos 2 9 + Sl, sin 2 9 
-sin(20){S x ,,S,,}/2), 



(135) 



-Jd 



cos(20) (s 2 x , cos 2 9 + S 2 , sin 2 9 - S 2 , 
+ sm(29){S x ,,S z ,}/2^ 
- sin(2</>) (cos 6{S X > , S y > } + sin 6{S y > , S z > }) 



(136) 



and 



H' d \<p?) = -^(2sin 2 0cos(20)[3m 2 - S ( S + l)]|^> 
+2 sin (9 {^m + a)Al m 



cr=±l 



x [cos 9 cos(2</>) + ia sin(20)] \ip" l+a ) 

+ F * m [t 1 + cos2 ) cos ( 2< ^) 

<7=±1 

+2iacos9sm(2(t>)]\ l fT +2a ] 



(138) 



where F* is defined by Eq. (17). We note that in this 
representation, both Hb and Hd preserve the global spin 
quantum number s, but allow Am = ±1, ±2 transitions. 
As a check on the algebra, we verified that for si = 1/2, 
Eq. (43) with J a = J c = is obtained from Eqs. (137) 
and (138). 



Local Hamiltonian 

With regard to the local spin anisotropy terms in the 
rotated coordinate system, we write 

H' a = -J a (Oi cos 2 9 + <D 2 sin 2 9 - ^% 3 )l39) 



and 



H' — -J, 



cos(2<£) (d sin 2 9 + 2 cos 2 9 

+ isin(20)O 3 -O 4 ) 
- sin(20)(o 5 cos + 6 sin , (140) 



where 



where {A, B} = AB + BA is the anticommutator. 



01 - £^„ (141) 

i=i 
2 

02 = (142) 

1=1 
2 

' Six' Six' Siz' 
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04 = J2 s l" ( 144 ) 

1=1 

2 

05 = ^^(-Six'Siy' + Si y 'Si x >), (145) 



and 



+ S iz ,S iy ,). (146) 

Local Hamiltonian matrix element components 

The operations of these interactions are given by 

+ E ^rVr +2CT ')), (147) 

<r'=±l 
i \ ^ ram |,„m+2<7\ 

<J=±1 

+ E ^t/i^w)), (148) 

ct,ct'=±1 

3 bD = J E fe>r +a > 

CT=±1 

- £ ^X^'l^)), (149) 

ct'=±1 

ct'=±i 

a=±l 

- E ^T/l^)), (150) 

cr,<7'=±l 

<j=±i 

+ E ^T; CT 'i^ + + 2 2 ')), (i5i) 

<r'=±l 

(7=±1 

- E ^X^Vr+^y (152) 



<r'=±l 



where 

M s m si = -4m 2 a s , Sl +4[s(a + l)-l](l-a a , ai ), 



(153) 



\rm,cr' _ r< -crcr m-(a +1) /2 n aa m+(a -1) /2 , r A \ 

Iy s, Sl - L 's+(cr' + l)/2,s 1 Ly s+(3(T'+l)/2,s 1 ^ 104 7 



P™f = 2A™(2m + a)a s , Sl , 
and 

ryx,n' fi— xa' — (a' + l)/2 jym+a 

n s,S! — L/ s+(cr' + l)/2,si s+(3ct'+1)/2,si 

-,-rECT'-(o-' + l)/2 



(155) 



+ °s+(3<x' + l)/2, Sl ^s+(<7' + l)/2, Sl i l 100 J 

where r? 5 , Sl , G™ t , AT*///, and L* Si are given by 

Eqs. (22) and (25)- (28), respectively. We note that for 
a' = ±l, N™>f =2H™>f. 



APPENDIX C 
Second order induction representation Hamiltonian 

In this appendix, we evaluate the corrections to the 
eigenstate energies second order in the four anisotropy 
interaction energies Jj for j = a, b, d, e. The operations 
of the rotated Hamilatonian TL 1 upon the eigenstates \ <p™) 
may be written as 

H'\<p?) = (E?^+E^)\^T)+ E W™s?\<P?+2*>) 

ct'=±1 

+ ^s>7 +a ) + Ks>7 +2a ) 

+ E ( x ™sr'\<p?X') 

2 2 CT ')), (157) 



<7,cr'=±l 



_i_ -ym.cr.cr' I m+2o 



where 



5, Si \ 5 



i(2m + a)Ar[sin(20)(^f 
-J d % Sl cos(2^ 



-ZiaJ^l 1 sin6>sin(2</>) 



(158) 



rricrm 



+^;e sl ( 1 + cos2 ) cos ( 2 ^ 

+2iaJ s / e 1 cos0sin(20)l, 



(159) 



1 



WJX (0, </>) = - [ J„ cos 2 + J e sin 2 9 cos(20)] 

+ iiV™f ' sin 2 0[J a - J e cos(20)], (160) 

= - ' (sin(20) [ J„ - J e cos(20)] 

-2iCTJ e sin0sin(20)Y (161) 



and 



--K" 



J a sin 



<T=±1 



+J e ((l + cos 2 6<) cos(2<; 



20 



sin(20))], (162) 



where E™s[°^ and are given by Eqs. (65) and 

(66), respectively, i?™ s 'f and ifg'£ are given by Eqs. 
(26) and (27), respectively, N™£* and -Rf;^' are given by 
Eqs. (154) and (156), respectively, and J^ 1 and Jj'^ 1 
are given by Eqs. (67) and (68), respectively. 

Second order eigenstate energies 

From Eq. (157), the second order eigenstate energies 
may be written as 



ET'W = — V a( 



cr=±l 



2 1 

+ 2 



V" 



+ E 



m,fr' |2 



^ J[2 + (2a + IK] 

S( 



p^m,(T,a' |2 



T> t7i 1 ^[2+(2.s + lH+a 7 i? 



J[2+ (2s + 1)ct']+ 2(775/ ■ 
For simplicity, we rewrite this as 



(163) 



™,(2) 
s,si 



S,Sl 



™,(2)V 

S,Sl 



E m,(2)U , r ,(2)V , £ .m,(2)W 

S,«l T S,Sl ^ S,Sl 

, F rn,(2)X , B m : (2)^ 

to sin 2 r , , i „ , 
^-g-[4s(s + l)-8m 2 -l] 

x(cos 2 0[J fc s ; Q Sl -cos(2^)4f] 2 
+ sin 2 (20)(J^) 2 ), 
[2s(s + l)-2m 2 -l] 



(164) 



(165) 



87B 



£ ,m,(2)W = 



E™,(2)* 

S,Sl 



x ([sin 2 + (1 + cos 2 0) cos(20)J d s ;: 1 ] 2 



4cos 2 0sin 2 (20)(J^ 1 ) 2 ), (166) 



16J 



j (j a - 3[J a cos 2 6 + J e sin 2 cos(20)]) , 

(167) 



/,7 ai (7£/J)sin 2 fl 
2J 

xf[J a - J e cos(2</>)] 2 cos 2 



+ sin 2 (20)J e ; 



(168) 



s,si 



9™ Sl {lB/J) 



( Jo sin 2 



64J 

+J e (l + cos 2 0)cos(20) 



+4cos 2 0sin 2 (2</>)J 2 



(169) 



where 



d m = - 



CM = 



( 5 2_ TO 2 )[(s _ 1) 2„ m 2 ] ^ i7? 2_ Mi 

{28-1) 

.2 2 
+ r ?s+2,si r ?s+l,si 

[( S +l) 2 -TO 2 ][(s + 2) 2 -TO 2 ] 

X (2s + 3) 

2 2 c 2 2N 

(2s- l) 2 -X 2 

x((2s-1)[(s-1)(s-2)+to 2 ] 

-to(2s - 3)a:) 

Vs+2, Sl Vs+i, Sl [(s + if - ™ 2 } 
(2s + 3) 2 - x 2 

x((2s + 3)[(s + 2)(s + 3) + to 2 ] 

-m(2s + 5)x), (170) 

+m 2 (6s 2 - 18s + 11) 
+s(s-l)(s-2)(s-3)] 

-4tox(2s - 3)(s 2 - 3s + 1 + m 2 )) (171) 



((2s + 3) [to 4 



(2s + 3) 2 - 4x 2 
+m 2 (6s 2 + 29s + 34) 
+ (s + l)(s + 2)(s + 3)(s + 4)] 

-4to:z;(2s + 5)(s 2 + 5s + 5 + to 2 )). (172) 

There is a remarkable amount of symmetry in the an- 
gular dependence of the eigenstate energies. We note that 



El 



) and E™si (0, <fi) have the same forms, dif- 
fering in the replacements of the interactions and 
Jj'* 1 with J a and J e , respectively, and with different 
overall constant functions. The same comparison can 



also be made with E' s "^ y (0,<j>) and ET,sT >v {0-, 4>)- In 
addition, we note that there is a remarkable similarity 



,(2)V, 



in the 0, 4> dependence of E, 



m,(2)W 



with that of the local 



spin anisotropy part of B 1 f t ^(6, 4>) in Eq. (73). 

APPENDIX D 

The contributions to the sth level crossing second order 
in the anisotropy interactions are calculated as indicated 
in Eq. (74), and are found to be 



( E s,(2) _ E s-h(2)\ 



B=-Js/i 



^a n (s,si)/ n (0,</>). 
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(173) 



Second order level crossing angular functions 



h{0A) 



sin 2 9 
J 



([J b -cos(2(/>)J d ] 2 cos 2 <9 

+J 2 sin 2 (2</>)), (174) 

h{6,4>) = ^([J 6 - J d cos(20)][J a - J e cos(20)] 

x cos 2 6» + J d J e sin 2 (20)]), (175) 

/ 3 (M) = ^([Ja-COs(20)J e ] 2 COS 2 

+J 2 sin 2 (20)), (176) 

1 / \ 2 

j (j 6 sin 2 9 + J d {\ + cos 2 9) cos(2</>)J 



/ 4 (6\0) - 



/ 5 (f,0) 



h(0, 



fr(0, 



I 2 

+4-|-cos 2 6>sin 2 (20), 

j[J fc sin 2 9 + J d (l + cos 2 9) cos(20)] 
x [J a sin 2 9 + J e (l + cos 2 0) cos(20)] 
+4^ cos 2 0sin 2 (20), 

j ( J a sin 2 9 + J e (l + cos 2 9) cos(2</>) 

+4^ cos 2 sin 2 (20), 



(178) 



(179) 



1 / \ 2 

0) = -[J a -3[ J a cos 2 6> + J e sin 2 6>cos(2</>)] J . 

(180) 

Second order level crossing coefficients 



The coefficients are given by 
(8s - 9) 

ai(s,si) — 
a 2 (s,s 1 ) -- 



2s ' 

[Q2,o(g) + 4si(.si + l)a 2 ,i(s)] 
s(2s + l)(2s + 3) 
-3(9 + s-23s 2 + 4s 3 + 12s 4 ), 
-9 + 8s + 4s 2 , 



fl2,o(s) 

a 2,l(s) 

d3{s, si) = a" (s,si) + a*(s,si), 



(181) 
(182) 

(183) 
(184) 



a^(s,si) 



a^(s,si) = 



[3-3s-3s 2 +4s 1 (s 1 + l)] 2 
2(2s + 3) 2 
(s - 1)[3 + 3s - 3s 2 + 4si(si + l)] 2 
2s(2s + l) 2 

[s(s + 2) -48^ + 1)] 
2(s + l)(s + 3)(2s + l) 2 (2s + 3) 2 



(185) 



<i(s) 

a 4 (s, si) 
a 5 (s, si) 

a5,o (s) 
a 6 (s, si) 

a^(s,si) 



<o(g)+ 481(81 + 1)0^(3) 

(2s + 5)(3s + l) 
(s + l)(s + 3)(51 + 114s 
+209s 2 + 302s 3 + 164s 4 + 24s 5 ), 
129 + 318s + 395s 2 
+442s 3 + 300s 4 + 72s 5 , 
(4s -3) 
8s ' 

3[a 5 , (s)+4si(si + 1)] 



4s(2s + l)(2s + 3) ' 
3 + 3s -5s 2 -4s 3 , 
ag(s, si) + a^(s,si), 
[3-3s-3s 2 +4si(si + l)] 2 
8(2s- l)(2s + 3) 2 
(s - 1)[3 + 3s - 3s 2 + 4si(si + l)] 2 



(186) 

(187) 

(188) 

(189) 
(190) 



(177) a£(s,si) = 



8s(2s- 3)(2s+ l) 2 
-1 



(191) 



a 6,l( s ) = 
a 6,2( s ) = 

a 7 (s,si) = 



96(2s - 3) (2s - l)(2s + l) 2 (2s + 3) 2 

ELo<»W[^i(«i + l)]" 
X (2s + 5)(4s+ l)(4s + 3) ' 

s(-5184 + 20898s + 123003s 2 

+81669s 3 - 304665s 4 - 521611s 5 

-168482s 6 + 163132s 7 + 124888s 8 

+ 12008s 9 + 256s 10 + 128s 11 ), (192) 

5832 + 6939s + 855s 2 + 17694s 3 

+22242s 4 - 1548s 5 - 3464s 6 

+3728s 7 + 1920s 8 - 512s 9 - 256s 10 , (193) 

4536 + 5787s - 3111s 2 + 4782s 3 

+ 16972s 4 + 5144s 5 - 864s 6 

+256s 7 + 128s 8 , (194) 

[s(s + 2)-4si(si + l)] 
4(2s+ l) 3 (2s + 3) 3 (2s + 5) 

x (4si(si + 1)(-1 + 38s + 24s 2 + 60s 3 ) 

+ (s + 1)(3 + 67s + 94s 2 + 44s 3 + 8s 4 )) . 

(195) 



By expanding the solutions in the crystal representa- 
tion to second order in the Jj , we have explicitly checked 
these expressions for si = 1/2, s = 1, and for si = 1, 
s = 1,2. We note that for si = 1/2, only a\ and a 4 are 
non- vanishing. 
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